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Pressure Induced Ferromagnetism in Cubic Perovskite SrFeC>3 and BaFeO; 
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The spin order in cubic perovskite SrFe03 and BaFe03 under high pressure is studied by den- 
sity functional theory (DFT) calculation with local spin density approximation plus Hubbard U 
(LSDA+f/). At ambient pressure, A-type and G-type helical spin orders are almost degenerate in 
BaFe03 whose lattice constant is 3.97 A. When the lattice constant is reduced to 3.85 A which is 
same as the lattice constant of SrFeC>3 at ambient pressure, G-type helical spin order becomes stable, 
being consistent with SrFeOs. This is because superexchange interaction is enhanced as compared 
with double exchange interaction. Phase transition from helical spin state to ferromagnetic state 
in both SrFeOs and BaFeOs takes place if the lattice constant is further reduced to 3.70 A. This 
is because reduced local spin moment weakens the contribution from superexchange interaction. 
Our result agrees with recent experimental result of BaFeOs under high pressure. Additionally, our 
calculation predicts that half-metal BaFeOs at ambient pressure will become a good metal under 
high pressure. 

PACS numbers: PACS numbers: 75.30.-m, 75.30.Et, 75.50.Bb, 75.40.Mg 



Helical spin order in cubic perovskite AFeOs (A=Ca, 
Sr, Ba), where Fe 4+ is in a high spin configuration d 4 , has 
attracted lots of research interest for their potential appli- 
cation in spintronic devices^— All of them present helical 
spin order below 115 K, 134 K, and 111 K for A=Ca, Sr, 
and Ba, respectively.— ~— The Fe3d electrons in these ma- 
terials can be divided into two classes: conducting and 
localized electrons. Three localized electrons occupy t2 g 
orbitals and one electron occupies double degenerated e g 
orbitals. The interaction between conducting electron 
and localized electron is described by Hund coupling. 
In addition, charge-transfer energy A defined as the en- 
ergy cost to move an electron from oxygen 2p orbital to 
Fe3d orbitals shows a negative valued— implying that 
metallic conduction mainly occurs on oxygen band. The 
helical spin order can be understood from the competi- 
tion of double exchange (DE) and superexchange (SE) 
interactions, the former and the latter of which favors 
ferromagnetism (FM) and anti- ferromagnetism (AFM), 
respectively 1 ^ 12 

It is well-known that electronic structure is tunable 
under high pressure. Phase transition from helical spin 
state to FM state in SrFeC>3 under the pressure of 7 GPa 
has been reported^ Very recently, the evolution of spin 
order in BaFeOs under pressure has been studied^ At 
ambient pressure, BaFeOs shows helical spin order, which 
changes to FM under very weak external magnetic field, 
~0.3 Ti£ The helical spin order is stabilized with increas- 
ing pressure, but FM finally becomes stable under pres- 
sure above 30 GPa. There is no structural phase transi- 
tion, since the cubic symmetry of BaFeOs preserves up 
to 50 GPa. The electrical resistance decreases under high 
pressure. 

Hydrostatic pressure P reduces the lattice constant a 
in AFe03. The compression of a leads to the increase 



of the hopping integral pda representing the hybridiza- 
tion between 02p and Fe3d orbitals. Since the DE and 
SE energies are roughly proportional to pda (ref. 11) and 
(pda) 4 , respectively, increasing pda gives the enhance- 
ment of SE as compared with DE. This explains the sta- 
bilization of helical spin order under pressure in BaFeOs, 
in which the G- and A-type helical spin orders are almost 
degenerate at ambient pressured However, this cannot 
explain the stabilization of FM under further pressure. 
Therefore, it is significant to study the electronic struc- 
ture of SrFeOs and BaFeOs under high pressure, which 
is helpful to understand the FM phase transition and the 
behavior of electrical resistance. 

In this paper, we perform density functional theory 
(DFT) calculations with local spin density approxima- 
tion plus Hubbard U (LSDA+J7) for both SrFe0 3 and 
BaFeOs. At ambient pressure, the lattice constant of 
SrFe0 3 and BaFe0 3 is 3.85 A and 3.97 A, respectively, 
and both materials show helical spin order. We find that 
BaFeOs represents stable G-type helical spin order that 
is the same order as SrFeOs, if the lattice constant is 
reduced to a = 3.85 A. We also find that FM is stable 
in both SrFeOs and BaFeOs if a is further reduced to 
3.70 A. Our result for BaFeOs agrees with recent exper- 
imental results under high pressured According to the 
DE model including oxygen orbitals and SE^i the in- 
crease of pda stabilizes FM when SE is unchanged. The 
SE energy is given by the product of neighboring spin 
moment with coupling parameter Jse- With decreasing 
a, Jse increases according to a relation Jse °c (pda) 4 , 
while the product of neighboring spin moment decreases 
because of the suppression of spin moment on Fe site. As 
a result, the SE energy may become less sensitive to a un- 
der high pressure. This view can explain the stabilization 
of FM obtained by our first-principles calculations. 
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FIG. 1. (Color online) The </> dependence of the total energy 
difference per unit cell, AE(cj>) = E(4>) — E(<f> = 0), obtained 
by LSDA+C/ for various lattice constant a shown as numbers 
in SrFeOs. (a) G-type helical spin order; (b) A-type helical 
spin order. At ambient pressure, a=3.85 A. For the definition 
of 6, see text. 
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FIG. 2. (Color online) The cj> dependence of the total energy 
difference per unit cell, AE(<f>) = E(<j>) — E(<f> — 0), obtained 
by LSDA+C/ for various lattice constant a shown as numbers 
in BaFeOa. (a) G-type helical spin order; (b) A-type helical 
spin order. At ambient pressure, a=3.97 A. For the definition 
of 6, see text. 



Non-collinear spin polarized calculations are performed 
by VASR±£ within LSDA+C/— in a primitive cell with 
10x10x10 k-point grid. For the non-collinear spin order 
calculation, the wave function is in the form of spinor and 
a generalized Bloch boundary condition is adopted, 1 -. 
The projector augmented wave pseudopotentials with 
Ceperly- Alder exchange-correlation with 500 eV energy 
cutoff is used in our calculation. The energy resolution 
is set to be 0.01 meV per unit cell. For band dispersions 
and density of states (DOS), 40x40x40 k-point grid is 
adopted for their self-consistent calculations. 

It is necessary to add on-site Coulomb interaction U 
and exchange interaction J to LSDA in order to obtain 
a helical spin order In the present study, we take a set 
of parameters, U = 4.0 eV and J = 0.9 eV, and for sim- 
plicity, we ignore the a dependence of U and J. These 
parameters are slightly larger than those used in a pre- 
vious paper (U = 3.0 eV and J = 0.6 eV)J^ However, 
similar to the previous study, this set also gives a propa- 
gating vector of helical spin order that is consistent with 
the observed one in SrFe03 at ambient pressure. Further- 
more, for BaFeOa a calculated lattice constant where the 



FM transition occurs agrees well with an observed one as 
will be shown below. We emphasize that the physics dis- 
cussed in the present paper is independent of the choice 
of the parameters. 

The propagating vector of A-type and G-type heli- 
cal spin orders is defined as q = —(<p, 0,0) and q = 
^ ((/>, <f>, <fj) , respectively. To find out the optimal value 
of <j), we define the (^-dependent total energy measured 
with respect to the energy of FM state, AE(<p) = 
E{(j>) - E((f> = 0). The energy in SrFe0 3 with differ- 
ent lattice constant is shown in Fig. [T] At ambient pres- 
sure (a=3.85 A), the G-type order with (f>=0.09 shows 
an energy minimum. This cf> is close to an observed 
one (0=O.112)i£ With decreasing a, the G-type order 
becomes more stable, but with further decreasing a a 
transition from helical spin state to FM state occures ir- 
respective of the type of order, i.e., the </> at minimum 
energy becomes zero when a changes from 3.75 A to 
3.70 A in both the A- and G-type orders. This result 
is qualitatively consistent with the experimental results 
of a transition from AFM to FM induced by pressured. 
An observed a at the transition pressure (7 GPa) is about 
a=3.78 A 13 ' 15 , which is slightly larger than the calculated 
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FIG. 3. (Color online) The spin moment M in BaFeOa with 
G-type helical spin order obtained by LSDA+[7 as a function 
of <f> for several values of a. 

one. The difference may partly come from the choice of 
U and J. We note that smaller U and J gives a larger a 
at the transition. 

The AE(4>) in BaFe0 3 with A-typc and G-type helical 
spin orders is shown in Fig. HJa) and Fig. HJb), respec- 
tively. At ambient pressure (et=3.97 A), the A-type and 
G-type orders are almost degenerate with small propa- 
gating vector With increasing pressure, i.e., reducing 
a from 3.97 A at ambient pressure to a=3.85 A, both the 
A-typc and G-type orders become more stable with larger 
4>. We note that at a=3.85 A the G-type order is more 
stable than the A-type order. The value of at mini- 
mal AE(4>) with a=3.85 A is almost the same as that in 
SrFeOs with the same a. As the case of SrFeOa, BaFe03 
shows FM if a is further reduced (see Figs.[5{a) and (b)). 
The pressure induced FM is thus common among ^4Fe03. 

The energy minimum of G-type helical spin order in 
BaFe03 (see Fig.[2fb)) shifts to the right-hand side when 
we reduce a from 3.97 A to 3.75 A. The shift of the en- 
ergy minimum can be explained by the competition be- 
tween DE and SE interactions. The reduction of a leads 
to the increase of pda. Since the DE and SE energies 
are roughly proportional to pda (ref. 11) and (pdc) 4 , re- 
spectively, the increase of pda gives the enhancement of 
SE interaction, comparing with DE. The enhanced SE 
stabilizes helical spin order at large 0. 

Figure [3] shows the calculated spin moment M as a 
function of for several values of a in BaFc03 with G- 
typc helical spin order. The magnitude of M decreases 
with decreasing a, and m , defined as the value of where 
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FIG. 4. (Color online) Density of states (DOS) of BaFe0 3 
in the FM state calculated by LSDA+17 with U = 4.0 eV 
and J — 0.9 eV. The positive side of the DOS denotes the 
up-spin DOS, while the negative side denotes the down-spin 
DOS. Dotted vertical line at zero energy represents the Fermi 
level. Black solid line and red solid lines represent the DOS 
of Fet2 g and Fee 9 orbitals, respectively. The partial DOS of 
02p orbitals is shown by solid blue line, (a) a— 3.97 A(ambient 
pressure), (b) a= 3.70 A(under pressure). 



M shows maximum, shifts toward larger 0. We find that 
the shift of m is consistent with the shift of the minimal 
AE(cj)) in Fig. H^b) when 3.97 A < a <3.75 A. This is a 
reasonable behavior, since it is expected that the energy 
is gained if the magnetization is large. 

With further reducing a to 3.70 A in Fig. [2][a), the 
global energy minimum shifts to a FM state. This means 
that DE finally comes over SE. We still can see the trace 
of a local minimum around 0=0.18. Correspondingly, the 
maximum of M is located around 0=0.18. It is interest- 
ing to notice that the transition to the FM state seems 
to be discontinuous, i.e., a first-order transition. This is 
clear from the E(<f>) curve for a=3.75 A, where there ap- 
pear two minima at 0=0 and 0.16. If the DE interaction 
is only between nearest-neighboring sites, E(4>) has only a 
single minimum. Therefore, the discontinuous transition 
indicates the presence of longer-range DE interactions. 
The change from helical to FM state between a=3.75 A 
and a=3.70 A obtained by our calculation agrees with the 
emergence of FM state in BaFe03 at a=3.73 A observed 
by experiment^. 
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FIG. 5. (Color online) Band structure of ferromagnetic 
BaFeC>3 calculated by LSDA+U. Solid and dotted lines repre- 
sent spin-up and spin-down bands, respectively, (a) a=3.97 A 
and (b) a=3.70 A. 



Now, a question is why the SE will fail to stabilize heli- 
cal spin order if the lattice constant reduced to a=3.70 A. 
The SE energy between nearest-neighbor localized spins 
is given by E SE = J SE J2<i,j> ( s * ' S j)' where J SE is the 
SE coupling parameter, Si is the spin operator at site i, 
< i,j > runs over the nearest-neighbor pairs, and (• • • ) 
represents an average of an static quantity. The SE en- 
ergy is, thus, proportional to the product of neighbor- 
ing spin moments. Since M decreases with decreasing 
a as shown in Fig. [3J (Sj • Sj) is also expected to de- 
crease. This will compensate the effect of the increase of 
Jse oc (pda) 4 , and thus the SE energy will become less a- 
dependent. As a consequence, the effect of DE overcomes 
the effect of SE, leading to the first-order transition to 
FM around a=3.70 A. 

To understand the mechanism of the reduction of M 
under high pressure, it is convenient to study DOS in 
BaFeC>3. Here, we focus on the cases with a=3.97 A and 
3.70 A, both of which show the energy minimum near 
or at cf>=0. Comparing DOS at a~3.97 A with that at 
a= 3.70 A shown in Fig. Eta) and Fig. 0Jb), we clearly 



see the increase of the energy width of DOS under pres- 
sure. We note that the increase of band width usually 
leads to smaller M through the suppression of DOS in 
the paramagnetic phase. 

The increase of band width also induces the change of 
DOS at the Fermi level. At a=3.97 A, only electrons with 
up spin contribute to the conductivity, and the DOS of 
down spin is almost zero at the Fermi level as shown in 
Fig. Ufa). This means that ferromagnetic BaFe03 is half- 
metallic at ambient pressure H At a=3.70 A, DOS at the 
Fermi level becomes nonzero and small part of t2 S DOS 
is occupied by down-spin electrons. This is consistent 
with the decrease of M. The band dispersions of BaFeOs 
in the FM state are shown in Fig. [5j Comparing the 
dispersions for a=3.97 A and a=3.70 A, we find that 
main changes under high pressure take place at the Y 
and X points. At the Y point, several spin-down bands 
descend and cross the Fermi level. At X point, one spin- 
up band lifts and one spin-down band descends. The 
downward shift of the spin-down bands is consistent with 
the emergence of extra t2 g electrons near the Fermi level 
in the down-spin DOS shown in Fig.[4fb) . In other words, 
a half metallic nature in BaFcOa at ambient pressure is 
lost under high pressure and becomes a good metal. This 
explains the observed suppression of electric resistance 
under high pressure.— 

In summary, the spin order of SrFe03 and BaFe03 
under high pressure is studied by DFT calculation with 
LSDA+J7. A transition from helical spin state to FM 
state in SrFe03 and BaFeOs under high pressure is re- 
produced. The transition is of the first-order type. The 
mechanism of the transition is associated with the evo- 
lution of DE energy and SE energy. The DE energy in- 
creases with decreasing a because of the enhancement of 
hopping integral pda. This favors FM. However, the gain 
of the SE energy competes with the DE energy, and as a 
result, the helical spin order is stabilized in the range of 
a=3.97~3.75 A. With further reducing a, the effect of SE 
fades out because of reduced local spin moment on Fe. 
Our calculated results agree with the recent experimen- 
tal result on BaFeOs. From the a dependence of DOS in 
the FM state, we predict that a half metallic character of 
ferromagnetic BaFeOs at ambient pressure is lost under 
high pressure. 
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